Surface dark solitons in nonlocal nonlinear media 
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We predict the existence of surface dark solitons at the interface between a self-defocusing nonlocal 
nonlinear medium and a linear medium. The fundamental and higher-order surface dark solitons 
can exist when the linear refractive index of the self-defocusing media is much larger than that of 
the linear media. The fundamental solitons are stable and the stabilities of higher-order solitons 
depend on both nonlocality degree and propagation constant. 
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Optical surface waves, a special type of waves that is 
localized at the interface between two different media, 
are widely used in sensing physical, chemical and biolog- 
ical agents to explore intrinsic and extrinsic properties of 
material surface. In the presence of nonlinearity, some 
kinds of surface solitons have been found theoretically 
and experimentally, such as in local Kerr median], 
waveguide arrays [3,IU, photorefractive media d,[^7and 
metamaterials [7]. Surface dark solitons existing at the 
interface between the local nonlinear media and the lin- 
ear media have been studied[l, Q. 

Nonlinear surface solitons can exist not only at the in- 
terface in local nonlinear media but also in nonlocal non- 
linear media. Varies types of surface solitons in nonlocal 
nonlinear media have been found and studied, such as 
bright surface fundamental solitons jljj , incoherent sur- 
face solitons surface dipoles [H,ll3| and surface vor- 
tices [l2|. In self-defocusing nonlocal media, the bright 
surface solitons and ring surface solitons have been pre- 
dicted and studied (HI, [l5j]. However, can surface dark 
solitons exist in self-defocusing nonlocal media? There is 
no study till now, though dark solitons in nonlocal self- 
defocusing medium were studied theoretically fl6l - fl8j and 
experimentallv fl9l - [22j . 

In this Letter, we address the existence of surface dark 
solitons (SDSs) at the interface between a self-defocusing 
thermal medium and a linear medium. We find that a 
SDS in a self-defocusing nonlocal media is identical with 
the half part of a dark soliton in the same bulk media 
when the linear refractive index of the self-defocusing 
media is much larger than that of the linear media. The 
similar relation between bright surface solitons and bulk 
solitons in nonlocal nonlinear media is firstly given in 
Ref.[23( under the same assumption. 

We start our analysis by considering a TE polarized 
laser beam with an complex envelope q propagating along 
the z axis in the vicinity of the interface formed by 
a self-defocusing thermal nonlinear medium and a lin- 



ear medium. The propagation of laser waves is de- 
scribed by the dimensionless (1+1)-D nonlocal nonlinear 
Schrodinger equation: 
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Here the transverse x and longitudinal z coordinates are 
scaled in terms of beam width and diffraction length, 
respectively. An denotes the nonlinear perturbation of 
refractive index in the self-defocusing thermal medium. 
rid > describes the difference of unperturbed refractive 
index between the self-defocusing thermal medium and 
the linear medium. The parameter w m stands for the 
characteristic length of the nonlocal material response. 
We use w m /wo represents the degree of the nonlocality, 
where wq is the beam width. 

For the TE polarized wave, the continuity condition 
for transverse fields is q(x = +0) = q(x = — 0). Follow- 
ing the experimental instance [10j, the continuity con- 
dition for the nonlinear refractive index is assumed as 
dAn/dx\ x= o — 0, which can be achieved when the in- 
terface between nonlinear and linear media is thermally 
insulating. We assume rid >> L which is easily satisfied 
in an actual physical system|23]. When rid >> 1 almost 
all the fields of laser beams will locates in the higher- 
index media [Toj - Tisl . l23j . so the boundary condition in the 
lower-index linear medium is q(x — > +oo) = 0. 

The asymptotic behaviors of SDSs when x — > — oo are 
similar to that of bulk dark solitons [l8j]. If the station- 
ary solution is in the form q(x,z) = W{x)e l ^ z , we have 
W(x — > —oo) — y/—/3. Therefore the propagation con- 
stant (3 must be negative. Consequently, one also have 
An(x — > — oo) = j3. In an actual thermal nonlinear sys- 
tem, An is proportional to the change of temperature 



2 




FIG. 1: (Color online) Profiles of fundamental SDSs for (a) 
w m = 1,3,5, /3 2 = l;(b)to m = 3, /3 2 = 1,1.5,2. Propagation 



of fundamental SDSs with (c)w„ 



i, F 



l;(d)«>„ 



ft — 1.5. Purple dash-dotted lines show the interface. 



due to the absorption of the wave energy. For a uni- 
form distribution of optical intensity I = \q\ 2 = |W| 2 (i.e. 
in the vicinity of x — > — oo), the distribution of tem- 
perature could not be uniform because the temperature 
gradient is necessary for the transfer of heat energy. As 
a result, the exact uniform background for surface and 
bulk dark solitons in thermal nonlinear media can not 
be achieved. A much broad beam is used to approx- 
imate the uniform background (same in local case) in 
experiments [19|, l20(. For the theocratical study of SDSs 
in this Letter, d&n/dx^^-oo = and dq/dx^^-oo = 
are used in numerical simulations for the convenience. 







w m =3 




w m =5 






FIG. 2: (Color online) (a)Profiles of second-order SDS with 
w m = 3,5, /3 2 — l.(b)The perturbation growth rate versus 
propagation constant for w m = 3, 5. Propagation of second 
order SDS for (c)w m = 3, f3 2 = 1; (d)w m = 3, /3 2 = 4.5; (e) 
w m = 5, /3 2 = 2.1; (f) w m = 5, /3 2 = 1. 



We search for stationary solutions for Eqs.([T]) and ([2]) 
numerically in the form q(x, z) = W(x)e t " z , where W(x) 
is a real function and the propagation Constance (3 < 0. 
The profiles of fundamental and higher-order SDSs are 
shown in FigsQ]|4j We find the fields of SDSs resides 
all in the self-defocusing media with higher refraction 
index. The profiles of SDSs are proved to be identical 
with the half part of corresponding bulk dark solitons, 
which are obtained numerically from Eqs.(fT]) for a bulk 
mcdium(— oo < x < oo). 

To elucidate the linear stability of SDSs, we search 
for perturbed solutions in the form q(x, z) — \W{x) + 
a(x, z)]e ll3z and An = N{x) + b(x,z), where \a(x, z)\ <C 
|W(x)| and |6(ir, z)\ <C |-ZV"(ac) ] are real perturbations. 
Substituting the perturbed solutions into Eqs.(fT]) and ©, 
one can get the linearized equations around stationary so- 
lutions W(x) and N(x) 
(i) in noninear media, x < 0, 



da 1 d 2 a 
l dz + 2dx^- Pa 
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(ii) in linear medis, x > 
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The perturbations a{x, z) and b{x, z) can grow with a 
complex rate a upon propagation. Following the method 
in Ref.[24j, one can get a linear eigenvalue problem for 
a. The eigenvalue problem has been solved numerically 
and the instability growth rate a r (real part of a) are 
shown in FigsJ2] -0J In addition, to confirm the results 
of linear stability analysis, we perform numerical simu- 
lations base on Eqs.([T]) and ([2]). The solutions of SDSs 
obtained numerical are used as the incident profiles with 
noise in the form q(x,z = 0) = T / F(a;)[l + p{x)], where 
p(x) is a random function with a Gaussian distribution 
and variance u 2 oise — 0.01. 

Figure Q] shows results for the fundamental SDSs with 
different nonlocality degrees and different propagation 
constants. As the nonlocality degree increases, the shoul- 
der of SDSs increases and moves away from the interface. 
The width of the intensity valley increases too. With in- 
creasing of P, the intensity of dark background increase 
and the width of valley decreases. The linear stability 
analysis shows the fundamental SDSs are always stable. 
The numerical simulations with noise- added incident pro- 
files, shown in Figsfllc) and (d), confirm the stability od 
SDSs. 

Figure [5] shows results for the second-order SDSs with 
different nonlocality degrees. There are two valleys in the 
second-order SDSs as shown in FigJ5Ja). As the nonlo- 
cality degree increases, the shoulder of the second-order 
SDSs increases and moves away from the interface. The 
width of two valleys increases together. The results of 
the linear stability analysis are shown in FigfJ^b). One 
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FIG. 3: (Color online) (a) Profiles of third-order SDS with 
w m = 2,5, f3 2 = 1. (b)The perturbation growth rate versus 
propagation constant for w m — 2, 5. Propagation of third- 
order SDS for (c) w m = 2, f3 2 = 2;(d) w m = 2, f3 2 = 4.6; (e) 
Wm = 5, P 2 = 1; (f) w m = 5, /3 2 = 2. 



can see that the stability of the second-order SDSs de- 
pends on both w m and /J. The numerical simulations 
with noise-added incident profiles are shown in FigsJ^c) 
- (f). In FigsHfc) and (e), the propagation are stable over 
1000 times diffraction length, as predicated by the sta- 
bility analysis[see Fig|2jb)]. In Figslljd) and (f) which 
are predicated unstable, the waves begin to oscillate over 
100 times diffraction length.. 

Figure [3] shows results for the third-order SDSs with 
different nonlocality degrees. There are three valleys in 
the third-order SDSs as shown in FigGJa). As the non- 
locality degree increases, the shoulder of the third-order 
SDSs increases and moves away from the interface. The 
width of three valleys increases together. The results of 
the linear stability analysis are shown in Fig|3fb). One 
can see that the stability of the third-order SDSs depends 
on both w m and f3. Comparing with the fundamental and 
second-order SDSs, the stable region for third-order are 
small. The numerical simulations with noise-added inci- 
dent profiles are shown in Figsj3jc) - (f). In Figsj3jc) and 
(e), the propagation are stable over 1000 times diffrac- 
tion length, as predicated by the stability analysis [see 
FigJ3fb)]. In FigsfJ^d) and (f) which are predicated un- 
stable, the waves begin to oscillate over 100 times diffrac- 
tion length.. 

Figure 2] shows results for the forth-order SDSs with 
different nonlocality degrees. There are four valleys in the 
forth-order SDSs as shown in Figj4ja). As the nonlocality 
degree increases, the shoulder of the forth-order SDSs 



increases and moves away from the interface. The width 
of three valleys increases together. The results of the 
linear stability analysis are shown in FigUJb). One can 
see that the stability of the forth-order SDSs depends on 
both w m and /?. Comparing with the fundamental and 
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FIG. 4: (Color online) (a)Profiles of forth-order SDSs for 
w m = 2,5, /3 2 — 1. (b)The perturbation growth rate ver- 
sus propagation constant for w m ~ 2, 5. Propagation of the 
forth-order SDSs for (c)w„ 
(e) w m = 5, /3 2 = 1; (f) w T . 



= 2, P 2 = 1; (d) Wr, 
= 5, P 2 = 3.4. 



2, p* = 5; 



second-order SDSs, the stable region for forth-order are 
much small. We have noticed that there are some 
noise in the curves in Figj4](b), the correctness 
of this results are still under examination. The 

numerical simulations with noise- added incident profiles 
are shown in FigsHJc) - (f). In Figsf^c) and (e), the 
propagation are stable over 1000 times diffraction length. 
In FigsHJd) and (f) which are predicated unstable, the 
waves begin to oscillate over 100 times diffraction length.. 

In summary, we introduced surface image dark solitons 
supported by the interface between self-defocusing ther- 
mal media and linear media. Such solitons integrate the 
unique features that are typical for surface waves sup- 
ported by nonlocal nonlinear interfaces with those ex- 
hibited by dark solitons existing in nonlocal nonlinear 
media. Results motivate further study to nonlocal dark 
solitons and surface solitons. 
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